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5 [i] Using recently developed techniques we estimate snow and sea ice thickness 

6 distributions for the Arctic basin through the combination of freeboard data from the 

7 Ice, Cloud, and land Elevation Satellite (ICESat) and a snow depth model. These data 

8 are used with meteorological data and a thermodynamic sea ice model to calculate 

9 ocean-atmosphere heat exchange and ice volume production during the 2003-2008 fall am 

10 winter seasons. The calculated heat fluxes and ice growth rates are in agreement with 

11 previous observations over multiyear ice. In this study, we calculate heat fluxes and ice 

12 growth rates for the full distribution of ice thicknesses covering the Arctic basin and 

13 determine the impact of ice thickness change on the calculated values. Thinning of the sea 

14 ice is observed which greatly increases the 2005-2007 fall period ocean-atmosphere heat 

15 fluxes compared to those observed in 2003. Although there was also a decline in sea 

16 ice thickness for the winter periods, the winter time heat flux was found to be less impacted by 

17 the observed changes in ice thickness. A large increase in the net Arctic ocean-atmosphere 

18 heat output is also observed in the fall periods due to changes in the areal coverage of 

19 sea ice. The anomalously low sea ice coverage in 2007 led to a net ocean-atmosphere heat 

20 output approximately 3 times greater than was observed in previous years and suggests that 

21 sea ice losses are now playing a role in increasing surface air temperatures in the Arctic. 
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25 1 . Introduction 

26 [ 2 ] Recent observations have shown a decline in Arctic 

27 sea ice areal coverage, freeboard, thickness, and volume 

28 [e.g., Stroeve et al, 2008; Farrell et al., 2000; Rothrock 

29 et al., 2008; Giles et al., 2008; Kwok et al., 2009] along 

30 with widespread environmental and climatic changes in the 

31 Arctic [Arctic Climate Impact Assessment, 2005]. These 

32 changes to the sea ice system have the potential to impact 

33 the Arctic climate by altering the radiation and heat budgets 

34 of the ocean and atmosphere. The degree to which the cold 

35 Arctic atmosphere is insulated from the relatively warm ocean 

36 is affected by the presence of a sea ice cover; the ocean- 

37 atmosphere heat flux can vary by nearly 2 orders of magnitude 

38 between open water and an ocean covered with thick sea ice 

39 for winter time conditions [Maykut, 1978]. This insulating 
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effect of sea ice makes the Arctic much colder than is typical of 40 
a maritime environment. The exchange of heat between the 41 
ocean and the atmosphere is also responsible for the growth 42 
of sea ice as heat lost from the ocean to the atmosphere is 43 
balanced by ice production. With thinner ice comes more 44 
heat exchange and faster ice growth which could potentially 45 
slow or reverse the observed losses in ice thickness. 46 

[ 3 ] The loss of sea ice may play a role in Arctic ampli- 47 
fication, wherein the Arctic region is expected to see a much 48 
greater share of warming as worldwide temperatures increase 49 
[. Manabe and Stouffer, 1980]. Modeling studies show that 50 
decreases in sea ice thickness and its areal coverage lead to 51 
increased ocean-atmosphere heat transfer. Due to the strong 52 
stratification of the Arctic atmosphere this heat is trapped 53 
near the surface leading to increased surface air temperatures 54 
[Boe et al., 2009]. In addition to modeling studies, observa- 55 
tions from buoy data have suggested that thinning of the sea 56 
ice cover during the 1979-1998 time period led to increases 57 
in surface air temperature through an increase in the ocean- 58 
atmosphere heat flux [Rigor et al., 2002], There remains, 59 
however, much uncertainty into how large a role recent 60 
changes in the sea ice cover have, and will continue to play, 61 
with regard to Arctic warming. Using reanalysis data, Serreze 62 
et al. [2009] found that losses in sea ice areal coverage have 63 
played a role in autumn surface air temperature increases in 64 
the Arctic. They also found that a winter warming signal may 65 
be beginning to emerge which they hypothesize may be due 66 
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Table 1. 

Input Parameters Used 

in This Study and Their Sources 

1 1 .2 

Symbol 

Description 

Source 

tl.3 

T a 

2 m air temperature 

ECMWF 

1 1 .4 

T d 2 

m dew point temperature 

ECMWF 

tl.5 

Po 

surface pressure 

ECMWF 

1 1 .6 

u 

10m wind speed 

ECMWF 

tl.7 

Cl 

cloud fraction 

MODIS 

tl.8 

T w 

sea surface temperature 

AMSR-E 

tl.9 

K 

snow depth 

snow model 

tl.10 

h f 

freeboard 

ICESat 

tl.l 1 

K 

ice thickness 

ICESat freeboard with snow model 


67 to delays in autumn freezeup and decreased ice extent and 

68 thickness in the winter. However, a major limitation in studies 

69 such as these has been the lack of a high-resolution, basin- 

70 wide sea ice thickness observational data set with which to 

71 adequately study the impact of sea ice thickness changes on 

72 the Arctic energy budget. 

73 [ 4 ] Recent satellite altimetry missions have provided the 

74 capability of obtaining basin-wide Arctic sea ice thickness 

75 measurements. In this paper, we use laser altimetry data 

76 from NASA’s Ice, Cloud, and land Elevation Satellite 

77 (ICESat) to estimate sea ice freeboard across the Arctic basin. 

78 The freeboard data are then combined with a snow depth 

79 model to estimate sea ice and snow thickness values for the 

80 Arctic at the high spatial resolution needed for studying the 

81 impact of sea ice on the energy budget. The sea ice thickness 

82 data are used with meteorological data and a thermodynamic 

83 sea ice model to study the impact of sea ice thickness changes 

84 on the ocean-atmosphere heat flux and ice growth rate over 

85 the 2003-2008 time period when significant changes to the 

86 Arctic sea ice cover took place. 

87 [ 5 ] The meteorological forcings, as well as the data sets 

88 and methodologies used to derive the sea ice thickness and 

89 snow depth are described in section 2. Sectioi/^describes 

90 the thermodynamic model used for determining the heat 

91 transfer through the ocean-ice-atmosphere system and cal- 

92 dilating the ice growth rate. The calculated heat fluxes, ice 

93 growth rates, and uncertainties are presented in section 4 and 

94 compared to results from previous studies. The role of 

95 observed thinning of the ice and snow covers in increasing the 

96 ocean-atmosphere heat flux is also discussed. Section 5 

97 expands the analysis to the full Arctic Ocean including 

98 nonice-covered regions. Section 6 summarizes the main 

99 conclusions of our study. 

100 2. Data Sets 

101 [6] In this section, we provide a description of the data 

102 sets and methods used to derive snow depth, sea ice thickness, 

103 and the meteorological parameters used in our analysis. These 

104 data sets are used in the following section to calculate the 

105 ocean-atmosphere heat flux and ice growth rate. No single 

106 sensor provides the requisite data, thus a combination of 

107 observation, model, and assimilated data is used. Table 1 

108 provides a summary of the input data sets with detailed 

109 descriptions provided below. Error estimates for each data 

110 set, along with the propagation of these errors into the 

111 calculated heat flux and ice growth rate, are addressed in 

112 section 5. We restrict our data set to the Arctic Ocean region 

113 shown in the shaded region of Figure 1 to avoid mixing 

114 high- and low-latitude sea ice regions in the analysis. 


2.1. Meteorological Data 115 

[ 7 ] Reanalysis data from the European Center for Medium- 116 

Range Weather Forecasts (ECMWF) ERA-Interim data set 117 
are used to provide the 2 m air temperature, 2 m dew point 118 
temperature, 10 m wind speed, surface pressure, and snow- 119 
fall. ERA-Interim combines observational and model data 120 
into an assimilated data set using the 4D-VAR method. Data 121 
is provided at 6 h time intervals with a spatial resolution of 122 
1.5° latitude by 1.5° longitude. 123 

[8] Cloud fraction is taken from the daily Moderate 124 

Resolution Imaging Spectroradiometer (MODIS) 1° x 1° 125 
global gridded product. A correction factor of 0. 1 has been 126 
added to all cloud fraction data to account for a bias in the 127 
Arctic region of the data set [ Ackerman et al., 2008]. Cloud 128 
fractions from MODIS, rather than ECMWF are used because 129 
of the anomalously high values found in the ECMWF data 130 
for this time period; the ECMWF cloud fractions were 131 
found to be approximately 30-40% higher than those from 132 
previously published observations [e.g., Lindsay, 1998]. 133 

[ 9 ] Sea surface temperatures are classified as the tem- 134 
perature of the top layer of water approximately 1 millimeter 135 
thick. They are taken from the daily 0.25° by 0.25° gridded 136 
product derived from ten-channel Advanced Microwave 137 
Scanning Radiometer-Earth Observing System (AMSR-E) 138 
brightness temperature data [Wentz and Meissner, 2004]. 139 
These sea surface temperatures are provided for ice-free 140 
areas to within 75 km of coastlines. The estimated error in 141 
the sea surface temperatures is 0.58 K [Wentz and Meissner, 142 

143 

2.2. Snow Model 144 

0 ] Snow depth on sea ice is modeled using a domain 145 




efmed by the 25 km AMSR-E grid. Snow depth on the 146 
model grid is determined by 147 

% = -V ■{ V -S)+aAF, 

dt p w 

where S is the average snow thickness in a grid cell 148 
(including both open water and ice covered areas), V is the 149 



Figure 1. Map of the region used in the analysis. The 
shaded region is defined as the Arctic Ocean in this study. 
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t2.1 Table 2. Time Periods Used in This Analysis Based on the 
t2.2 Availability of ICESat Data 3 


t2.3 

Campaign Name 

Period 

Days of Operation 

t2.4 

ono3 

Oct 1 to Nov 18 2003 

49 

t2.5 

ON03 1 

Oct 1 to Nov 8 2003 

39 

t2.6 

ON03 2 

Oct 15 to Nov 18 2003 

35 

t2.7 

FM04 

Feb 17 to Mar 21 2004 

34 

t2.8 

ON04 

Oct 3 to Nov 8 2004 

37 

t2.9 

FM05 

Feb 17 to Mar 24 2005 

36 

t2.10 

ON05 

Oct 21 to Nov 24 2005 

35 

t2. 1 1 

FM06 

Feb 22 to Mar 27 2006 

34 

t2. 12 

ON06 

Oct 25 to Nov 27 2006 

34 

t2.13 

MA07 

Mar 12 to Apr 14 2007 

34 

t2. 14 

ON07 

Oct 2 to Nov 5 2007 

37 

t2. 15 

FM08 

Feb 17 to Mar 21 2008 

34 


t2. 16 a The ON03 campaign has been subdivided into two campaigns, ON031 
t2.17 and ON03_2, for better temporal comparison with other fall ICESat 
t2.18 campaigns. 


150 ice velocity vector, a, is the ice concentration, p s is the snow 

151 density, p w is the density of water, and F is the snowfall (in 

152 snow water equivalent). The snow depth is initialized each 

153 year on 15 September before the summer minimum sea ice 

154 extent, the initial snow cover on multiyear ice and the snow 

155 density values are taken from the climatology of Warren et al. 

156 [1999]. The daily AMSR-E sea ice concentrations at each 

157 grid point are specified at the start of each day and remain 

158 constant throughout the day. Daily snowfall at each model 

159 grid point is estimated using the liquid water equivalent from 

160 the ECMWF ERA-Interim reanalysis data similar to the 

161 method used by Kwok and Cunningham [2008], Ice velocity 

162 for each grid point is determined from AMSR-E 89 GHz data 

163 using the wavelet analysis algorithm of Liu and Cavalieri 

164 [1998], The model is run each year during the fall through 

165 spring periods to estimate the snow depth over the time 

166 period covering each ICESat measurement campaign. 

167 2.3. ICESat Data 

168 [n] ICESat measures the surface elevation using a 

169 1064 nm laser altimeter [Zwallv et al., 2002], Spatial cover- 

170 age of the Arctic Ocean is provided up to 86°N with a 170 m 

171 shot-to-shot spacing and a footprint size of approximately 

172 70 m. The cloud filtering parameters described by Kwok 

173 et al. [2007] are first used to filter out low-quality data 

174 which has been affected by atmospheric forward scattering. 

175 The elevation data from ICESat are used to determine the sea 

176 ice freeboard, hf, which is here defined as the height of the 

177 snow and ice layer above the local sea surface. Freeboard 

178 data is collected only in areas where the ice concentration 

179 determined from AMSR-E is greater than 30%. The ICESat 

180 data products are of Release 428, which include orbit and 

181 attitude determination as well as detector saturation correc- 

182 tions for the time periods studied here. Freeboard is found 

183 from the ICESat elevation data through the use of sea surface 

184 tie points following the method of Kwok et al. [2007], 

185 [ 12 ] Due to the approximately 70 m footprint size of 

186 ICESat, some sea surface tie points used in the retrieval of 

187 freeboard from ICESat data are expected to be biased due to 

188 contamination of snow and ice within the footprint. Corn- 

189 parisons of ICESat data with coincident high-resolution air- 

190 borne laser altimetry data have shown this can be problematic 

191 with a freeboard bias of up to 9 cm observed in one study 


[Kurtz et al., 2008]. Corrections to account for biases due to 192 
snow and ice within sea surface tie point footprints have been 193 
proposed by Kwok and Cunningham [2008] and Kwok et al. 194 
[2009] and are applied here in the determination of free- 195 
board. The correction for snow depth biases are taken from 196 
Kwok and Cunningham [2008] which relates the albedo 197 
dependence of snow depth to the surface reflectivity mea- 198 
sured by ICESat. An additional correction to account for 199 
remaining residual biases due to contamination of snow and 200 
ice within the ICESat footprint is taken from Kwok et al. 201 
[2009]. 202 

[ 13 ] The temporal sampling of ICESat is limited to the 203 

times shown in Table 2 which restricts our analysis to time 204 
periods when ICESat data is available. Throughout we will 205 
refer to ICESat campaigns by their campaign name shown in 206 
Table 2, the first two letters of the campaign name refer to 207 
the months of measurement while the numerals refer to the 208 
year (e.g., ON03 for the October-November 2003 campaign). 209 
The length of the ON03 campaign made it suitable to split 210 
into two subcampaigns for the purposes of comparing the 211 
heat flux and ice growth rates between years. The ON03_1 212 
campaign is at a similar time of year to the ON04 and ON07 213 
campaigns while the ON03_2 campaign is at a similar time 214 
of year to the ON05 and ON06 campaigns. The FM04, FM05, 215 
FM06, and FM08 ICESat campaigns occurred during roughly 216 
the same time of year while the MA07 campaign occurred 217 
later in the ice growth season than all other campaigns. 218 

2.4. Sea Ice Thickness and Snow Depth 219 

[ 14 ] The sea ice thickness, h h is calculated by assuming 220 

local hydrostatic balance and is given by 221 


where hf is the height of the snow and ice layers above 222 
the water level, h s is the snow depth, p w = 1024 kg m -3 223 
is the density of sea water, p t is the density of sea ice 224 
taken to be 915 kg m 3 [Weeks and Lee, 1958; Wadhams 225 
et al., 1992], and p s is the density of snow. p s is taken to 226 
be changing with time following the climatological values 227 
compiled by Warren et al. [1999], it varies from a min- 228 
imum of 260 kg m 3 in early October to a maximum of 229 
330 kg m at the end of the winter ICESat campaigns. 230 

[ 15 ] The large difference between the spatial resolutions 231 
of the freeboard (approximately 70 m) and snow depth 232 
(25 km) data sets leads to ambiguities when combining 233 
these data to estimate sea ice thickness. Due to the nonlinear 234 
dependence of the heat flux values on snow and ice thickness 235 
(an example of which can be seen in Figure 2 for typical 236 
winter time conditions), it is necessary to use a high spatial 237 
resolution estimate of the thickness values to properly include 238 
the contributions of thin, young ice regions which can be 239 
present in any area due to ice dynamics. Kurtz et al. [2009] 240 
found that the mean heat flux and ice growth values ealeu- 241 
lated for the Arctic basin using the full 70 m spatial reso- 242 
lution of ICESat were approximately one-third higher than 243 
those calculated using 25 km mean thickness values. There- 244 
fore, the method developed by Kurtz et al. [2009] for com- 245 
bining low-resolution snow depth data with high-resolution 246 
freeboard data is used to estimate the snow and ice thickness 247 
distributions for each of 25 x 25 km grid cells in the Arctic 248 
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Sea ice thickness [m] 

Figure 2. Plot of the dependence of the ocean-atmosphere 
heat flux on sea ice thickness for snow-free and snow- 
covered sea ice using typical winter time conditions in the 
Arctic. Input parameters are as follows: air temperature of 
-25°C, cloud fraction of 0.5, wind speed of 6 m/s, relative 
humidity of 0.9, and no shortwave flux. 


249 containing available ICESat freeboard data. The method is 

250 based on an observed linear relationship between freeboard 

251 and snow depth for thin ice. The linear relationship between 

252 freeboard and snow depth applies to points with a freeboard 

253 less than a certain cutoff value, fb cuto ff. fb cuto ff is defined as 

fbcutoff = 0.69 (h s ) + 0.22 (h f ) + 5.10, 

254 where (/;.,) is the mean snow depth of the region which is 

255 given by the 25 km resolution snow depth model, • /(,) i> the 

256 mean freeboard of the ICESat data line within the 25 km 

257 snow depth grid cell, and the units of the constant, 5.10, are 

258 in cm. A constant snow depth is used for thick ice (where 

259 hf > fb cutoff ) and is given by 



hsthick = 1-03 {hi) + 0. 


260 where the units of.the constant value, 0.83, are also in cm. 

261 h s is thus given by 


h \ hSihick 


n 1 hf < fbcutoff 

'l \TDcutoffJ 

( htfyhick hf ft) cuto ff 


The discrete ICESat ice and snow thickness data points are 
assumed to represent the thickness distribution in each 
model grid cell, and the heat flux and ice growth values are 
calculated for each individual ice thickness data point in a 
grid cell containing a valid number of measurements. Heat 
transfer between the ocean, ice, snow, and atmosphere is 
governed by the temperature of each system, the tempera- 
tures of the ocean and atmosphere are specified, while the 
temperature profiles of the ice and snow are calculated. The 
temperature of the ocean layer in contact with the ice is 
taken to be near the freezing point of seawater at T b = 
271.35 K, while the surface air temperature and other rele- 
vant meteorological parameters are taken from the ECMWF, 
AMSR-E, and MODIS data discussed in section 2. Tem- 
perature gradients are mainly vertical, therefore disregarding 
horizontal heat fluxes the temperai 
snow and ice layers is governed 
heat diffusion equations 
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10 J J kg~' K _1 and c ice = 2.1 x 291 
1 0 ’ J kg -1 K 1 are the specific heats of ice and snow, and 292 
k s = 0.31 W m _1 KT 1 and k, = 2.04 W nT 1 K _1 are the 293 
thermal conductivities of snow and sea ice, respectively, which 294 
are empirical values obtained from Maykut and Untersteiner 295 
k[1969]. A more recent study by Sturm et al. [2002] also 
found the effective thermal conductivity for snow to be 
approximately 0.3 W m 1 K 1 . The numerical scheme used 
to solve equations 2 and 3 follows the three-layer model of 
Semtner [1976] with parameterizations for the individual 
heat flux terms described in detail below. 

[n] The resultant mean surface air temperature, ocean- 
atmosphere heat flux, and ice growth rates used in sections 4 
and 5 are the model average values over each ICESat mea- 
surement time period. They were calculated by running the 
thermodynamic model with 6 h time steps over each specific 
time period shown in Table 2. The initial temperature pro- 
files of the snow and ice layers were determined by first 
setting the system in thermodynamic equilibrium then run- 
ning the model over a one week time period prior to the start 
of each campaign shown in Table 2. 


296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 


3.1. Heat Flux Parameterizations 


312 


262 Here hf is taken from the ICESat data set, and /?, is then 

263 calculated for each freeboard data point using equation 1. 

264 The ice thickness distribution for each 25 x 25 km grid cell 

265 is then estimated from the approximately 70 m resolution 

266 ice thickness data. A minimum of 70 freeboard points (about 

267 half the grid cell coverage) are required for the determination 

268 of the ice thickness distribution in each grid cell. 

269 3. Thermodynamic Sea Ice Model 

270 [ 16 ] The ocean-atmosphere heat fluxes and ice growth 

271 rates are calculated here through the use of a thermodynamic 

272 model with inputs from the data sets described in section 2. 


[is] The various heat flux terms are calculated by solving 313 
the energy balance equation to find the surface temperature, 314 
T 0 , based on the method of Maykut [1978]. The energy 315 
balance equation at the surface is 316 

F r + F b — F E + F s + F e + F c = 0, (4) 

where F r is the net absorbed surface shortwave flux, F L the 317 
incoming longwave flux, F E the emitted longwave flux, F s 318 
the sensible heat flux, F e the latent heat flux, and F c the 319 
conductive heat flux. A positive flux is defined as being 320 
toward the surface while a negative flux is away from the 321 
surface. 322 
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323 [ 19 ] The net absorbed shortwave flux, F r , can be written 

as 


F,. = F ro (l-a)(l-*o), 


(5) 


324 where F ro is the shortwave flux reaching the surface, a is 

325 the surface albedo, and i 0 is the percentage of shortwave 

326 radiation which passes through the surface and into the water. 

327 For snow covered ice a is 0.8 and io is 0. For ice with a 

328 negligible snow cover (<1 cm thick is treated here as snow 

329 free) a is a function of ice thickness, h t , and calculated using 

330 the empirical relation between ice thickness and albedo 

331 described by Weller [1972]. i 0 is estimated from radiative 

332 transfer calculations described by Mavkut [1982], 

333 [ 20 ] Many parameterizations of the F r and F L radiative 

334 flux terms have been proposed in the literature. Key et al. 

335 [1996] analyzed various schemes and found that the short- 

336 wave parameterization scheme of Shine [1984] and the 

337 downwelling longwave parameterization scheme of Mavkut 

338 and Church [1973] perform well for Arctic conditions. F,. o is 

339 calculated here following Parkinson and Washington [1979] 

340 by applying the cloudiness factor of Laevastu [1960] to the 

341 empirical equation of F r for clear skies described by Shine 

342 [1984]. The downwelling longwave parameterization scheme 

343 of Maykut and Church [1973] is used to calculate F L . 

344 [ 21 ] The emitted longwave radiation, F E , is given by 

F e = eaTl (6) 


345 where e is the longwave emissivity of the surface layer taken 

346 to be 0.99, a is the Stefan-Boltzmann constant, and T 0 is tin 

347 temperature of the surface layer. 

348 [ 22 ] The turbulent fluxes are calculated using bulk aero 

349 dynamic formulas following Pease [1987] 


F s = pc p C s u(T a - T 0 ) 


F e = pLC e u(q a - q 0 


\Zj 

r> 


(V) 


( 8 ) 


so the number of grid points is reduced as needed to maintain 373 
computational stability. For the case of ice with a thickness 374 
less than 22 cm, the “zero layer” method of Semtner [1976] is 375 
used to determine the vertical temperature profile, the snow 376 
and ice layers are treated as a single system that maintains 377 
thermodynamic equilibrium with the external conditions at 378 
all times. 379 

[ 23 ] The ocean-atmosphere heat flux is defined as the net 380 
heat transferred from the ocean to the atmosphere, or ~F C . 381 
For open water areas, the individual heat flux terms are 382 
calculated using the above relations for F,., F L , F E , F s , and 383 
F e with suitable changes to a, i 0 , T 0 , and L. The surface 384 
albedo of open water is taken to be 0.08 while i 0 is the 385 
amount of shortwave energy passing through the ocean 386 
mixed layer which is calculated to be 0.2 based on the results 387 
of Maykut and Perovich [1987] for a 30 m mixed ocean 388 
layer. The latent heat of sublimation, L, is replaced by the 389 
latent heat of vaporization which is 2.5 x 10 6 J kg -1 . The 390 
surface temperature, T 0 , is replaced by the ocean surface 391 
temperature, T w . T w is taken to be constant at 271.35 K for 392 
ice-covered regions. The net ocean-atmosphere heat flux is 393 



- F7- F l - F s - F e . 


(9) 


ic Ice Growth Rate 


3.2. Therm 

[ 24 ] Ablation and accretion of ice at the bottom of the sea 
ice layer occurs when there is an imbalance between the 
conductive flux through the bottom of the ice ( F cn ) and the 
rergy from the water to the ice ( Fq ). The thermo- 
ic basal ice growth rate is calculated as 



394 

395 

396 

397 

398 

399 


;(*■«- 4 )> 


dhj 1 . 

17-eP-- F ' 


( 10 ) 


400 


350 where p is the air density, c p = 1004 J is the specific 

351 heat of air at constant pressure, C s 1 0 3 and C e = 2 x 

352 10 3 are the sensible and latent heat transfer coefficients, 

353 respectively, for neutrally stratified air and are adjusted for 

354 unstable conditions following Hack et al. [1993], u is the 

355 average wind speedjS = liS3 x 10 6 J kg~' is the latent heat of 

356 sublimation, and q is the specific humidity. The conductive 

357 flux, F c , is calculated by following the three-layer model of 

358 Semtner [1976]. Three vertical grid points are used: one in the 

359 snow layer, and two evenly spaced grid points in the ice layer. 

360 The surface energy balance equation (equation 4) can now be 

361 rewritten through substitution of the parameterizations for F r , 

362 F l , F e , F„ F e , and F c . The surface temperature-dependent 

363 terms in the surface energy balance equation are linearized to 

364 determine the temperature change of the surface layer for each 

365 time step. A time step of 6 h is used to coincide with the 

366 temporal resolution of the input ECMWF meteorological data 

367 described in section 2. Due to the coarse resolution of the 

368 temperature grid, a forward differencing scheme is used to 

369 calculate the conductive fluxes across the snow and ice layers 

370 and find the temperature profile, which is assumed to be linear 

371 between interior grid points. The forward differencing scheme 

372 is stable for vertical grid points with h t > 22 cm and h s > 14 cm, 


where Q t = 3.02 x 10 J m is the heat of fusion of ice, 

Fq is estimated to be 2 ± 1 W nT 2 from the results of 401 
Steele and Boyd [1998], and F cn is the conductive flux 402 
through the lowest ice grid point. The thermodynamic growth 403 
rate is calculated only to estimate the mean rate of ice growth 404 
for the observed ICESat thickness distributions, it is not 405 
used to change the thickness of the ice with time. 406 

4. Results for the Ice-Covered Arctic Ocean 407 

[ 25 ] The results presented in this section are for the sea ice 408 

covered region of the Arctic Ocean containing valid ICESat 409 
data. The ocean-atmosphere heat fluxes and ice growth rates 410 
represent approximately a monthly mean value for the study 411 
region. 412 

4.1. Heat Flux and Ice Growth in Regions Containing 413 
ICESat Data 414 

[ 26 ] Changes in the percentage distribution of different 415 

Arctic sea ice thickness classes over the 2003-2008 time 416 
period are shown in Figure 3 for both the fall and winter 417 
time periods. A general thinning of the ice cover is observed 418 
due to the loss of ice with thickness greater than 3 m. This is 419 
consistent with recent studies showing much of the older, 420 
thicker multiyear ice cover of the Arctic being replaced with 421 
thinner first year ice [Maslanik et al., 2007; Comiso et al., 422 
2008], Using similar data sets and methods, Kwok et al. 423 
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Arctic fall 


Arctic winter 



2003 2004 2005 2006 2007 2008 
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Arctic fall 




Figure 3. Distribution of ice thickness classes over the Arctic basin for the (a) fall and (b) winter ICESat 
campaigns, (c and d) The mean effective insulation of the snow plus sea ice cover in terms of an equiv- 
alent thickness of snow-free sea ice is also shown. The dark colored bars in Figures 3c and 3d represent 
the sea ice contribution, while the lighter colored bars represent the snow depth contribution. 


424 [2009] showed a comparable thinning of the Arctic sea ice 

425 cover with an overall decrease in the mean thickness over 

426 the same time period. The sea ice thickness results shown 

427 here differ from those of Kwok et al. [2009] due mainly to 

428 differences in the sea ice density used ( Kwok et al. [2009] used 

429 pi = 925 kg nT 3 while this study uses p t = 915 kg nT 3 ). 

430 Wadhams et al. [1992] summarize the results of numerous 

431 field measurements from the 1950s through the 1970s which 

432 suggest the mean density of sea ice is typically within the 


range 9 1 0-920 kg nT 3 for first year ice and 9 1 0-9 1 5 kg m -3 433 
for multiyear ice. However, whether the density of sea ice has 434 
changed with time due to changing ice conditions is an 435 
important, but unknown factor in the determination of sea ice 436 
thickness. Errors in the calculated heat flux and ice growth 437 
rates due to uncertainty in sea ice density are discussed in 438 
section 5. Figure 3 also shows the changes that occurred to 439 
the mean effective insulation of the sea ice cover over this 440 
time period. The effective insulation is defined here as the 441 
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442 thermal insulating strength of the snow plus sea ice layer in 

443 terms of an equivalent thickness of snow-free sea ice, it is 

444 calculated as h e jf = h t + jh s . The effective insulation of the 

445 fall ice pack decreased significantly in 2005 then remained 

446 relatively constant. The loss in the effective insulation 

447 during the fall periods is associated mainly with thinning of 

448 the sea ice rather than a loss of snow. During the winter 

449 time periods, the effective insulation stayed relatively con- 

450 stant until 2008 when it decreased by approximately 1 m 

451 (Figure 3). This decrease in the winter of 2008 is due to 

452 thinning of both the sea ice and snow covers which is 

453 associated with the large loss in multiyear ice and record 

454 minimum sea ice extent observed in 2007. 

455 [ 27 ] The percentage of ice within a given ice thickness 

456 class and the area weighted heat flux values for the various 

457 thickness classes are shown in Table 3 . Also shown in Table 3 

458 are the following mean input parameters: 2 m air temperature, 

459 cloud fraction, wind speed, and the calculated surface tem- 

460 perature. The calculated values are for areas where free- 

461 board data from ICESat were available which can be seen in 

462 Figures 4 and 5. Areas without ICESat data were not con- 

463 sidered in the analysis in this section. 

464 [ 28 ] Table 3 shows that over half of the ice production and 

465 ocean-atmosphere heat flux (~F C ) in the ice-covered regions 

466 of the Arctic Ocean occurred over areas with an ice thick- 

467 ness less than 80 cm. In particular, open water and newly 

468 reftozen leads with an ice thickness less than 10 cm accounted 

469 for nearly one-third of the ocean-atmosphere heat flux and ice 

470 production within ice-covered areas. The thickest ice (>1.6 m) 

471 is the dominant ice type and was found to make up 50-60% 

472 of the total observed ice in the Arctic. Yet, the thickest ice 

473 accounted for only 20-30% of the observed ice production 

474 and ocean-atmosphere heat flux. The basin wide averaged 

475 ice growth rate was generally higher in the winter than in 

476 the fall, this was due to the lower surface air temperatures 

477 and increased area of first year ice during the winter periods. 

478 The percentage contribution of each thickness class to ice 

479 production and heat flux varied due to the changing ice 

480 thickness distributions and input meteorological parameters. 

481 The net radiative flux showed the highest variability of the 

482 radiative, turbulent, and conductive heat fluxes. However, 

483 if we exclude the anomalous MA07 time period from corn- 

484 parison (which had a higher net radiative flux due to the 

485 increased shortwave flux of the later spring period) the net 

486 radiation was almost constant and varied by only 4 W m 2 . 

487 The loss of radiative energy by the atmosphere was observed 

488 to be much stronger over areas of thick ice rather than thin 

489 ice. The sensible heat flux was quite variable with variations 

490 of 8 W nY 2 seen during the study period. It acted to transfer 

491 heat from the surface to the atmosphere over relatively wann, 

492 thin ice (//,- < 0.4 m), while over ice thicker than 0.4 m, it 

493 transferred heat from the atmosphere to the surface. Overall, 

494 the sensible heat flux was positive owing to the large areas 

495 of thick ice in the Arctic, this resulted in a net sensible heat 

496 gain by the ice. The latent heat flux varied by 2 W nY 2 for 

497 all time periods and was generally a source of small but 

498 steady heat input to the atmosphere. 

499 [ 29 ] The input forcings and calculated heat flux values 

500 from this study are compared with results and observations 

501 from studies by Lindsay [1998], Maykut [1982], and Persson 

502 et al. [2002] in Table 4. The results shown in Table 4 for 


this study represent the mean over sea ice 2.75-3.25 m thick 503 
to best correspond with the observations conducted on 504 
multiyear ice floes in the comparison studies. The computed 505 
heat fluxes and forcing parameters derived in this study are 506 
within the range of observational values, with the exception 507 
of the sensible heat flux and surface air temperature, which 508 
were found to be slightly higher during the fall periods. We 509 
also compare our results for ice growth rates with those 510 
observed during the Surface Heat Budget of the Arctic 511 
Ocean (SHEBA) experiment. Perovich et al. [2003] studied 512 
basal ice growth rates for a 1.75 m thick multiyear ice floe 513 
(“Quebec site”) which grew to about 2.25 m thick between 514 
early October and March 1998. They report growth rates of 515 
0.10-0.30 cm d 1 in the fall and 0.25-0.50 cm d 1 in the 516 
winter (at comparable times to the fall and winter ICESat 517 
campaigns shown in Table 2). For a similar ice thick- 518 
ness class (ice of thickness between 1.75 and 2.25 m), 519 
we obtained similar Arctic-wide growth rates of 0.19- 520 
0.32 cm d 1 (mean 0.24 cm d *) in the fall and 0.27- 521 
0.44 cm d 1 (mean 0.33 cm d 1 ) in the winter. These 522 
comparisons demonstrate reasonable agreement between 523 
our derived results and observations from previous studies. 524 
The maj or advantage of the remote sensing data sets used here 525 
is that it is now possible to calculate the ocean-atmosphere 526 
heat flux and ice growth rate for all ice-covered areas of the 527 
Arctic. Table 3 thus expands on the knowledge from pre- 528 
vious observational studies by providing information over 529 
the full range of ice thickness classes of the Arctic Ocean. 530 

[ 30 ] Maps of the mean effective insulation, surface air 531 

temperature, ocean-atmosphere heat flux, and ice growth 532 
rate are shown in Figure 4 for the fall time periods and 533 
Figure 5 for the winter time periods. Figures 4 and 5 show 534 
that there was great spatial and temporal variability in the 535 
effective insulation, air temperature, heat flux, and ice growth 536 
rate during the study period. An analysis of the variability in 537 
the heat flux and ice growth rate, due to losses in the effective 538 
insulation coupled with changes in the meteorological for- 539 
cings, is the subject of section 4.2. 540 

4.2. Analysis of Heat Flux and Ice Growth Variability 541 

[ 31 ] The mean values for the ocean-atmosphere heat 542 

fluxes and ice growth rates in Table 3 do not show a clear 543 
correlation between an increased ocean-atmosphere heat 544 
flux/growth rate and the observed decrease in ice thickness 545 
and snow depth derived from the ICESat and snow model 546 
data sets. This follows since the observed heat flux also 547 
depends on the various meteorological forcings with the 548 
surface air temperature playing the largest role. Since sur- 549 
face air temperatures in the Arctic tend to be highly variable, 550 
it is likely that any trend in the heat flux values over this 551 
short 5 year time period is masked by the natural variability 552 
caused by variations in the surface air temperature. 553 

[ 32 ] The goal of this section is to better understand the 554 
causes of the variability that occurred over the study period. 555 
That is, we seek to determine whether the observed vari- 556 
ability of the heat flux and ice growth is due mainly to 557 
changes in meteorological conditions, changes in ice and 558 
snow thickness, or uncertainties in the input parameters. 559 
First, we first determine the uncertainty in the heat flux and 560 
ice growth rates through estimation of the errors in the input 561 
parameters. Next we run the thermodynamic model for each 562 
time period using constant meteorological forcings to focus 563 
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0.1 


Table 3. Thickness Distribution Averages, Ice Production, and Heat Flux Values Over the Ice-Covered Regions of the Arctic Ocean 3 


t3 .2 

Thickness Category 

ON03 1 

ON03J2 

FM04 

ON04 

FM05 

ON05 

t3 .3 




Percentage of Ice in 

Each Thickness CaU 

t3 .4 

0-0.1 m 

1.3 

1.3 

1.4 

1.3 

1.7 

1.5 

t3 .5 

0. 1-0.2 m 

0.6 

0.7 

0.9 

0.8 

1.2 

1.1 

t3 .6 

0.2-0. 4 m 

2.6 

2.8 

3.4 

2.9 

4.2 

4.7 

t3.7 

0.4-0. 8 m 

9.9 

11.6 

13.0 

10.7 

14.3 

17.5 

t3 .8 

0.8-1. 6 m 

21.9 

23.7 

28.3 

22.0 

26.0 

29.1 

t3 .9 

1. 6-3.0 m 

29.5 

28.9 

26.2 

28.8 

23.0 

26.3 

t3.10 

>3.0 

34.2 

31.0 

26.6 

33.6 

29.7 

19.8 

t3.11 








t3.12 




Net Radiation F r + F L — F E 

(W fM~ 2 ’ 

t3.13 

0-0.1 m 

-1.1 

-1.3 

-1.6 

-1.2 

-1.6 

-1.4 

t3. 14 

0. 1-0.2 m 

-0.3 

-0.3 

-0.4 

-0.3 

-0.5 

-0.5 

t3.15 

0.2-0. 4 m 

-0.9 

-1.1 

-1.3 

-1.0 

-1.5 

-1.8 

t3.16 

0.4-0. 8 m 

-3.1 

-3.9 

-4.2 

-3.3 

-4.6 

-5.9 

t3. 17 

0.8-1. 6 m 

-6.1 

-6.9 

-8.3 

-6.3 

-7.5 

-8.8 

t3.18 

1.6-00 

-14.3 

-13.7 

-12.8 

-14.6 

-12.1 

-11.3 

t3.19 

Total 

-25.8 

-27.1 

-28.6 

-26.7 

-27.7 

-29.8 

t3 .20 








t3 .2 1 





Sensible Heat Flux F s (W m 2 ) 

t3 .22 

0-0.1 m 

-2.7 

-3.0 

-4.8 

-3.0 

-4.3 

-3.4 

t3 .23 

0. 1-0.2 m 

-0.2 

-0.3 

-0.5 

-0.2 

-0.5 

-0.3 

t3 .24 

0.2-0. 4 m 

0.0 

-0.2 

-0.6 

-0.2 

-0.4 

-0.2 

t3 .25 

0.4-0. 8 m 

1.0 

0.8 

0.3 

0.8 

1.0 

1.5 

t3 .26 

0.8-1. 6 m 

2.9 

2.6 

3.2 

2.6 

3.3 

3.6 

t3 .27 

1.6-00 

7.0 

6.0 

6.5 

6.9 

6.3 

5.5 

t3 .28 

Total 

7.9 

6.0 

4.1 

6.9 

5.3 

6.6 

t3 .29 








t3 .30 





Latent Heat Flux F e (W 

m~ 2 ) 

t3 .3 1 

0-0.1 m 

-1.0 

-0.9 

-1.3 

-1.0 

-1.3 

-1.2 

t3 .32 

0. 1-0.2 m 

-0.1 

-0.1 

-0.1 

-0.1 

-0.2 

-0.1 


FM06 


ON06 


MA07 


ON07 


FM08 


0.33 

0.34 

0.35 

0.36 

0.37 

0.38 

0.39 

0.40 

0.41 

0.42 

0.43 

0.44 

0.45 

0.46 

0.47 

0.48 


0.2-0. 4 m 
0.4-0. 8 m 
0.8-1. 6 m 
I.6-00 
Total 


0-0.1 m 
0. 1-0.2 m 
0.2-0. 4 m 
0.4-0. 8 m 
0.8-1. 6 m 
I.6-00 
Total 


1.5 
1.1 

3.6 

15.2 

31.3 
23.7 
23.5 


1.4 
0.9 
3.7 

16.4 
28.6 
27.6 
21.3 


I. 3 
0.7 
2.7 

II. 7 

29.5 
28.4 

25.6 



Conductive Heat Flux F C (W m ~) 


22.2 


Ice Growth Rate (cm month ) 


1.5 

0.9 

3.8 

15.4 

29.8 

30.5 
18.1 



1.3 

0.7 

2.7 

14.2 

31.9 

32.1 

17.2 


-1.4 

-0.3 

- 1.0 

-4.8 

-9.7 

-12.7 

-29.9 


-3.8 

-0.4 

-0.4 

0.5 

3.4 

6.0 

5.3 


- 1.1 

- 0.1 

- 0.2 

-0.5 

- 0.6 

- 0.2 

-2.7 


5.2 

7.1 

6.0 

6.7 

6.2 

4.9 

4.7 

6.3 

0.6 

1.2 

1.0 

1.0 

0.9 

0.5 

0.5 

0.8 

1.3 

2.3 

2.2 

1.8 

1.9 

1.1 

1.3 

1.7 

2.8 

4.2 

4.9 

4.0 

4.9 

2.4 

3.1 

4.7 

3.8 

4.5 

5.4 

4.9 

5.6 

3.6 

4.1 

6.8 

7.3 

5.7 

5.6 

4.8 

5.8 

3.6 

4.7 

7.0 

21.0 

25.0 

25.2 

23.2 

25.3 

16.2 

18.4 

27.2 


t3.49 

0-0.1 m 

4.1 

im 4 .1 

V 6.5 

4.5 

6.1 

5.1 

5.8 

5.3 

4.2 

4.0 

5.4 

t3 .50 

0. 1-0.2 m 

0.4 Wr 

0.5 

0.9 

0.5 

1.0 

0.8 

0.8 

0.7 

0.4 

0.5 

0.6 

t3 .5 1 

0.2-0. 4 m 

0.9 

k 1 T 

1.8 

1.1 

1.9 

1.8 

1.5 

1.6 

0.9 

1.0 

1.4 

t3 .52 

0.4-0.8 m 


V^-5 

3.5 

2.1 

3.3 

3.8 

3.1 

3.8 

2.0 

2.2 

3.8 

t3 .53 

0.8-1. 6 m 

2.0 

7 2.8 

4.2 

2.3 

3.3 

3.8 

3.6 

4.0 

3.1 

2.5 

5.4 

t3 .54 

1.6-00 

2.5 

r 3.3 

4.6 

3.1 

3.2 

3.1 

3.2 

3.4 

3.8 

2.2 

5.6 

t3 .55 

Total 

Til. 6 

14.7 

21.6 

13.6 

18.8 

18.5 

18.0 

18.9 

14.4 

12.5 

22.2 

t3 .56 


V 











t3 .57 


y 



Mean Input Parameters 







t3 .58 

{ T a ) (K) 

253.8 

250.2 

244.5 

252.9 

248.3 

251.7 

249.1 

250.8 

253.3 

257.8 

246.6 

t3 .59 

( T s ) (K) 

251.8 

248.2 

242.7 

251.0 

246.3 

250.2 

247.0 

249.7 

251.8 

256.0 

245.1 

t3 .60 

(Cl) 

0.64 

0.58 

0.42 

0.61 

0.48 

0.58 

0.44 

0.67 

0.55 

0.65 

0.45 

t3 .6 1 

(u) (m s ) 

6.2 

5.5 

6.2 

6.0 

6.3 

6.2 

6.5 

6.2 

6.8 

6.7 

5.7 


t3 .62 “The heat fluxes and ice production rates for the different ice thickness categories have been weighted by the percentage of ice within each respective 
t3 .63 thickness category. 


564 exclusively on how the observed changes to the sea ice and 

565 snow thickness distributions affected the heat flux and 

566 growth rates across the Arctic ice pack. 

567 4.2.1. Sensitivity to Input Parameter Uncertainties 

568 [ 33 ] We now estimate the sensitivities and uncertainties 

569 in the heat flux and growth rate due to variations in the 

570 input parameters. To determine the impact of variability in 


the input parameters on the heat flux and ice growth rate, 571 
the thermodynamic model was run multiple times to simulate 572 
variations in each individual parameter separately over a 573 
range of values. The goal was to calculate the sensitivities of 574 
the heat flux (^) and ice growth rate ( dgr g^ th ) to the input 575 
parameters (x), and estimate an uncertainty value by multi- 576 
plying the sensitivity by the estimated uncertainty, a x . Sea- 577 
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Figure 4. Map of the effective insulation, snow depth, and air temperature parameters and the calculated 
ocean-atmosphere heat fluxes and ice growth rates for the fall measurement periods. 


VU 

l and used in the estimation 
; fluxes and ice growth rates in 
of the calculated sensitivities 
the fall and winter time 
the following discussion, 
ertainties are assumed to be from a 
ncess. All other error sources are not 
ke net error estimates a Fc and a growth 
are RSS errors calculated from the 


578 sonal sensitivities were calcula 

579 of the uncertainties of the 1 

580 section 4.2.2. Average vali 

581 and estimated uncertainties 

582 periods are shown 

583 only the freeboard 

584 zero mean random pr 

585 well constrained, thus 

586 presented in Table 5 

587 individual error terms. 

588 [ 34 ] Estimating uncertainties for the meteorological input 

589 parameters is challenging since errors in the ECMWF Interim 

590 surface air temperature, and wind speed for the Arctic have 

591 not been adequately determined at this time. For sea ice 

592 covered regions, the ECMWF meteorological parameters are 

593 modeled assuming a uniform snow-free 1.5 m thick ice slab, 

594 ice concentration is considered using a blend of model and 

595 observation data [Stark et al., 2007]. As shown in Figures 4 

596 and 5, the assumption of a unifonn effective ice thickness 

597 of 1 .5 m is typically not valid which may impact the ECMWF 

598 model results. The uncertainties in the ECMWF data depend 

599 not only on the model accuracy, but also on the quantity and 


quality of observations used in the assimilation which can 600 
vary considerably in time and space. Here we estimate the 601 
uncertainties in these values by assuming that they represent 602 
50% of the maximum observed variability of the areal mean 603 
across similar time periods. For example, the mean surface 604 
air temperature of the ice-covered Arctic, (T a ), varied from 605 
253.3-257.9 K between the ON03_1, ON04, and ON07 606 
campaigns leading to an observed variability of 4.6 K and 607 
an estimated uncertainty of 2.3 K. Similarly, uncertainties of 608 
0.6 m/s were estimated for the wind speed. Lupkes et al. 609 
[2010] compared ECMWF Interim near surface air tern- 610 
peratures and wind speeds to data from several ship cruises 611 
in the late summer in the Arctic and found a warm bias of 612 
1.5-2 K in the Interim temperature data set and near zero 613 
error in the wind speed. While this bias in the summer data 614 
may not apply to the fall and winter time periods used in 615 
this study, it suggests that our uncertainties for the surface 616 
air temperature and wind speed may be a reasonable estimate. 617 
However, the uncertainty in the surface air temperature may 618 
vary regionally as it depends on the number of observations 619 
used in the assimilation. Additionally, the low resolution of 620 
the ECMWF data could potentially lead to errors near the ice 621 
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Figure 5. Map of the effective insulation, snow depth, and air temperature parameters and the calculated 
ocean-atmosphere heat fluxes and ice growth rates for the winter and early spring measurement periods. 



ud fractions are estimated to 
d on a study by Ackerman 


622 edge. Errors in the MODIS c 

623 be 0. 1 for the Arctic region b; 

624 etal. [2008], 

625 [ 35 ] Errors in the ice thickness^tnd snow depth input 

626 parameters are due to uncertainties in the freeboard, snow 

627 depth, and density values. Errors in the freeboard were 


assumed to be unbiased (after the corrections for biases due 628 
to snow and ice contamination were applied) but estimated 629 
to have a random normally distributed error of <Jj h = 5 cm 630 
[ Kwok and Cunningham, 2008]. a p . is estimated to be 631 
10 kg/m 3 which represents the range of expected densities 632 
for sea ice between 0.3 and 3 m thick [Kovacs, 1996], cr is 633 


t4. 1 Table 4. Comparison of.Heat Flux and Forcing Parameters for the 
t4.2 Mean of All 2.75-3.25 m Thick Ice Areas With Observations' 1 


Table 5. Sensitivity of the Ocean- Atmosphere Heat Flux and Ice t5.1 
Growth Rate to Variations in the Input Parameters 3 t5.2 


t4.3 

Parameter 

This Study 
(2.75-3.25 m 
Ice Only) 

L98 

M82 

P02 

t4.4 

Net radiation 

-22 (-23) 

-24 (-26) 

-23 (-18) 

-20 (-20) 

t4.5 

F s 

12(11) 

8 (4) 

12(5) 

5(5) 

t4.6 

Fe 

0(1) 

1 (0) 

0 (-2) 

-1 (1) 

t4.7 

Fc 

11 (11) 

- 

11 (14) 

6(10) 

t4.8 

Ta( K) 

248 (252) 

241 (250) 

242 (249) 

251 (250) 

t4.9 

u { ms ; 

6(6) 

4(4) 

5 (5) 

5(7) 

t4. 10 

Cl 

0.5 (0.6) 

0.5 (0.6) 

- 

- 


t4. 1 1 “Observations are from Lindsay [1998] (L98), Maykut [1982] (M82), and 
t4.12 Persson et al. [2002] (P02). Values from M82 are taken from the 3 m ice 

t4. 1 3 thickness results. Values for the fall time periods are in parentheses, while 

t4.14 those for the winter are not. 


X 

cr x 

Heat Flux (W m 2 ) 

Growth Rate 
(cm month -1 ) 

t5.9 

dF c 

Ox 

„ 9Fc 

dx 

d growth 
dx 

dgrowth 
° x dx 

T«(K) 

2.3 

1. 1(1.0) 

2. 5(2.3) 

0.9(0.8) 

2.1(1.8) 

t5 .5 

Cl (%) 

10 

0.02(0.01) 

0.2(0. 1) 

0.02(0.01) 

0.2(04) 

t5 .6 

u (m s 1 ) 

0.6 

0.8(0. 8) 

0.5(0. 5) 

0.7(0.7) 

0.4(0. 4) 

t5.7 

ft>s! (cm) 

5 

0.3(0. 3) 

1 .6(1.5) 

0.3(0. 3) 

1.4(1. 3) 

t5.8 

h s (cm) 

5 

0.02(0.01) 

0.09(0.04) 

0.02(0.01) 

0.08(0.04) 

t5.9 

Pi (kg m 3 ) 

10 

0.1(0.1) 

0.9(0. 8) 

O.l(O.l) 

0.8(0. 7) 

t5 .1 0 

Ps (kg m 3 ) 

100 

0.01(0.01) 

0.7(1 .2) 

0.01(0.01) 

0. 6(1.0) 

t5.1 1 

Fb (W m 2 ) 

1 



0.9 

0. 9(0.9) 

t5.12 

°F C 



3.3(3.2) 



t5.13 

O growth 





2. 8(2. 7) 

t5 .1 4 


“Results for the winter time periods are in parentheses. t5. 1 5 
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634 estimated to be 100 kg/m 3 based on the variability of p s in 

635 the climatology of Warren et al. [1999]. Uncertainties and 

636 sensitivities due to variations in the density of sea water, 

637 dew point temperature (humidity), and surface air pressure 

638 are small and not considered here. Errors in the snow depth 

639 are unknown and estimated to be 5 cm here, but this value 

640 will be shown to be of small importance in the following 

641 discussion. 

642 [36] Table 5 shows that most of the uncertainty in both the 

643 heat flux and ice growth values is due to the relatively large 

644 uncertainty estimated for T a with lesser contributions due to 

645 uncertainty associated with sea ice freeboard, cloud fraction, 

646 wind speed, snow density, and ice density. Errors due to 

647 snow depth uncertainties are minor and contribute little to 

648 uncertainties in the heat flux and growth rates since errors in 

649 the snow depth are nearly canceled by the corresponding 

650 retrieval errors in ice thickness. Essentially, 1 cm of snow has 

651 an effective insulation of kj/k s = 6.5 cm of ice, while a 1 cm 

652 error in snow depth leads to a corresponding error of p p w Z P p ‘ ~ 

653 6.5 cm in ice thickness which makes errors due to snow depth 

654 uncertainties small. In this assessment, errors in the calculated 

655 mean heat flux and ice growth rate values for the Arctic are 

656 primarily due to errors in T a . However, changes in the cloud 

657 cover and associated incoming longwave radiation can also 

658 lead to changes in the surface air temperature which cannot 

659 be studied with a simple model such as the one used here. 

660 Aside from the impacts to surface air temperature, cloud 

661 cover changes are not a strong source of variability in the 

662 sensitivity of the ice growth rate and heat flux values. To 

663 better estimate the errors in the heat fluxes and ice growth 

664 rates calculated here, additional studies of the error in the 

665 ECMWF data for T a in the Arctic during the fall and winter 

666 time periods are needed. The next largest source of error is 

667 due to freeboard uncertainties, these errors are due to instru- 

668 mental uncertainties and set a lower limit for tha total 

669 uncertainty in the calculated heat flux and ice growth rate. 

670 4.2.2. Heat Flux Variability in Ice-Covered Regions 

671 [ 37 ] The sensitivity results for the various meteorological 

672 forcings shown in Table 5 demonstrate that changes in T a 

673 are much more dominant than Cl and u in affecting vari- 

674 ability in the calculated heat fluxes and ice growth rates. 

675 Variability in the surface air temperature is therefore one of 

676 the main factors that must be considered in analyzing the 

677 observed variability in the ocean-atmosphere heat flux and 

678 ice growth rate. Figure 6 shows the mean ocean-atmosphere 

679 heat flux and ice growth rate for the ice-covered Arctic 

680 Ocean over the different time periods as well as the corre- 

681 sponding mean surface air temperatures. The observed heat 

682 fluxes and growth rate values can be seen to primarily change 

683 with variations in the surface air temperature. However, the 

684 changes in the ON05, ON06, ON07, and FM08 time periods 

685 are disproportionate compared to earlier changes in T a . The 

686 ON05 and ON06 heat fluxes were much higher than those 

687 observed during the ON03_2 time period despite the higher 

688 surface air temperatures. Similarly, the winter FM08 time 

689 period has a higher heat flux than the FM04 time period 

690 despite a higher surface air temperature of 2.1 K. Figure 3 

691 shows that there was a significant change in ice thickness 

692 distribution and an associated large decline in the effective 

693 insulation during these time periods. The percentage of ice 

694 with a thickness greater than 3 m experienced the greatest 

695 decline beginning around the fall of 2005 and this was 


accompanied by an increase in the percentage of 0.4-1. 6 m 696 
ice in the fall and 0.8-1. 6 m ice in the winter. As shown in 697 
Figure 2, the ocean-atmosphere heat flux is sensitive to 698 
changes in the percentage of thin ice, especially for ice less 699 
than approximately 1 m thick. The percentage of the thin- 700 
nest ice classes (<0.4 m) did not change significantly over 701 
the 2003-2008 time period, however this value is reported 702 
for the ice-covered Arctic only and does not take into 703 
account the large changes in open water and loss of ice area 704 
for the entire Arctic also observed during this time period. 705 

[38] The FM05 and FM06 time periods have similar mean 706 
growth rates, heat fluxes, and surface air temperatures 707 
(Figures 6b and 6d) even though there was a decline in the 708 
percentage of thick ice during this time and a decline in 709 
mean ice thickness of 38 cm. The decrease in the percentage 710 
of the ice >3 m thick was compensated by an increase in the 711 
percentage of ice 1. 6-3.0 m thick (Figure 3b). Since the 712 
ocean-atmosphere heat flux and ice growth are much less 713 
sensitive to changes for ice in this thickness range it appears 714 
that variability in heat flux and ice growth during these 715 
winter time periods was dominated more by variability in 716 
the surface air temperature. The MA07 heat flux and growth 717 
rate is much lower than the other winter time periods, this is 718 
likely due to the higher surface air temperatures resulting 719 
from the later date of data collection as well as thicker ice 720 
cover due to the longer time available for sea ice growth. 721 

[ 39 ] The full effect of the observed increase in the ocean- 722 

atmosphere heat flux due to a thinning of the ice and snow 723 
cover is difficult to quantify since the ocean-atmosphere heat 724 
flux and surface air temperature are coupled. The ocean- 725 
atmosphere heat flux will increase with decreasing tempera- 726 
ture and vice versa until an equilibrium is reached between 727 
the surface heat flux and other factors (such as atmospheric 728 
energy transport) which determine the surface air tempera- 729 
ture. Nevertheless, to investigate the effect of changes in the 730 
snow and ice thickness distribution on the observed heat 731 
flux values (independent of changes due to meteorological 732 
conditions), we ran the thermodynamic model for the ice 733 
and snow thickness distributions for each individual time 734 
period using the same fixed meteorological conditions. 735 
Figure 7 shows the ocean-atmosphere heat flux differences 736 
for the individual time periods under the same meteorological 737 
conditions relative to the first campaign of the fall or winter 738 
season. This shows that thinning of the sea ice and snow 739 
covers led to potential ocean-atmosphere heat flux increases 740 
of nearly 6 W nT 2 for the fall 2005-2007 time periods 741 
compared to the 2003 time period (an increase of approxi- 742 
mately 40% over the heat flux observed in ON03_1). Despite 743 
the similarly large decrease in the effective insulation 744 
observed in ON05 and FM08 (Figure 3), the FM08 ocean- 745 
atmosphere heat flux would only be 2 W nT 2 higher than 746 
FM04 under equivalent meteorological conditions (an increase 747 
of approximately 10% from the observed heat flux in FM08), 748 
but this is also within the uncertainty of the values. 749 

[ 40 ] The results show that the observed thinning of sea ice 750 
during the 2005-2008 time period led to large increases in 751 
the ocean-atmosphere heat fluxes for the subsequent fall 752 
periods. The increased ocean-atmosphere heat flux likely 753 
impacted the surface air temperatures and may have played 754 
a part in the surface air temperature anomalies observed 755 
during this same period by Serreze et al. [2009]. The winter 756 
results suggest that despite losses in ice thickness and 757 
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Arctic fall 


Arctic winter 



Figure 6. The mean ocean-atmosphere heat flux, basal ice growth rate, and 2 m air temperature for ice- 
covered regions during the Arctic fall and winter seasons. 
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effective insulation, growth of the sea ice and the addition of 
snow over the fall and early winter limited increases to the 
winter heat flux. The MA07 results show a lower equivalent 
heat flux than FM04 which is due to the additional time for 
growth for the thin ice classes which reduces the overall 
heat flux. The FM08 results suggest that an increase in the 
ocean-atmosphere heat flux may be beginning to appear in 
the winter due to the large decrease in ice and snow thickness 
(effective insulation), however this cannot be fully deter- 
mined here due to uncertainties in the input parameters. 


4.3. East and West Arctic Differences 

[ 41 ] Sections 4.1 and 4.2.2 showed that ice thickness and 
energy exchange for the ice-covered regions of the Arctic 
Ocean experienced changes for the 2003-2008 time period, 
however certain regions of the Arctic were impacted dif- 
ferently than others. Here we discuss the regional impact of 
such changes by dividing the Arctic into two regions, East 
Arctic (0°-180° longitude) and West Arctic (180°-360° 
longitude), for the purpose of studying the regional vari- 
ability of ice thickness, energy exchange, and ice growth. 
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Arctic fall 



Arctic winter 



Figure 7. Ocean-atmosphere heat flux differences for the 
different time periods under the same meteorological condi- 
tions, differences are relative to the first campaign of the 
season. The error bars for the heat flux differences are taken 
from the combined uncertainties from the freeboard, snow 
depth, snow density, and ice density uncertainties discussed 
in section 4.2.1. 

778 The regional ice thickness distributions, mean surface air 

779 temperatures, and mean growth rates are discussed. The mean 

780 growth rate and heat flux tenns are used interchangeably 

781 here since the two are closely related. 

782 [ 42 ] Figure 8 shows a large decline in the amount of thick 

783 ice (>3 m) in both regions during the fall periods, with the 

784 East Arctic showing a particularly steep decline in 2005. 


Much of the ice of thickness greater than 3 m was replaced 785 
by ice 0.4-1. 6 m thick, with large increases in the 0.2-0. 8 m 786 
ice thickness class in the East Arctic. Both regions experi- 787 
enced similar variabilities in the surface air temperature, but 788 
differences in growth rate variabilities can be seen between 789 
the eastern and western Arctic regions due to differences in 790 
the ice thickness distribution. In 2005 and 2006 the East 791 
Arctic region experienced sharp increases in the ice growth 792 
rate/heat flux compared to the ON03_1 period (which had a 793 
lower surface air temperature) due largely to the increased 794 
amount of 0.2-0. 8 m thick ice. The West Arctic region 795 
experienced similar, but less prominent, increases in the ice 796 
growth/heat flux in 2005 and 2006 due to the loss of thick 797 
ice >3 m. 798 

[ 43 ] Figure 9 shows the regional thickness distributions, 799 

ice growth rate, and surface air temperature for the winter 800 
periods. The East Arctic winter time periods also experienced 801 
a general decline in the percentage of thick ice >3 m while 802 
the West Arctic did not see large changes in the ice thick- 803 
ness distribution until 2008. Despite losses in the thickest 804 
ice category as well as the overall mean ice thickness, the 805 
ice growth rate/heat flux is similar for the respective regions 806 
with similar surface air temperatures. Thus, as was observed 807 
in section 4.2 for the ice-covered Arctic, most of the winter 808 
time variability in ice growth rates appears to be due to 809 
changes in surface air temperature rather than due to changes 810 
in the ice thickness distribution. 811 

5. Results for the Full Arctic Ocean 812 

[ 44 ] Section 4 showed changes to the ocean-atmosphere 813 

heat flux and ice growth rate for areas containing ICESat 814 
data. We now extend the analysis to the full Arctic Ocean, 815 
including open water areas, to better place the results into 816 
context given the large changes in sea ice areal coverage 817 
over the time period. 818 

[45] In this section, the heat flux and ice growth rates are 819 

calculated for nonice-covered areas by using sea surface 820 
temperature data described in section 2. Areas with an ice 821 
concentration greater than 0 and less than 30% were treated 822 
initially as open water, but with a sea surface temperature at 823 
the freezing point of sea water. For the nonice-covered areas, 824 
the ice growth rate and ocean-atmosphere heat flux were 825 
calculated at 6 h time intervals. If the sea surface temperature 826 
was at the freezing point the ice was allowed to grow in 827 
thickness and the growth rate was approximated from the net 828 
surface heat flux and equation 10, if the sea surface temper- 829 
ature was greater than freezing point of sea water then the ice 830 
thickness and growth rates were set to 0. Without the insu- 831 
lation of a sea ice cover, the net surface heat flux tended to be 832 
much larger than that from the ice-covered regions. However, 833 
the rate of ice growth rate is not directly proportional to the 834 
net surface heat flux in nonice-covered areas because of the 835 
limitation that ice will only grow once the surface temper- 836 
ature has reached the freezing point. 837 

[46] To determine the net heat output and ice production of 838 
the Arctic Ocean, we first grid the heat flux and ice growth 839 
rate data onto a 25 km polar stereographic grid. Gaps in the 840 
gridded data were filled in through the use of a Gaussian 841 
smoother with a 20 km length scale (following Kwok et al. 842 
[2009]). Ice-covered and nonice-covered areas were filled 843 
in independently using their respective data sets. The pole 844 
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Figure 8. Fall time period ice thickness distributions, mean basal ice growth rates, and mean surface air 
temperatures for the ice-covered east and west Arctic regions. 


845 hole north of 86 degrees was not filled in due to the large 

846 uncertainty introduced in interpolating the data over such a 

847 large region. The total area of the Arctic Ocean considered 

848 in this section for all time periods is 6.47 x 10 6 km 2 . The 

849 net surface heating rate and net ice volume production are 

850 this area value multiplied by the ocean-atmosphere heat flux 

851 and ice growth rates, respectively. Results for the net surface 

852 heating rate and ice volume production as well as the areal 

853 coverage of ice and nonice areas are shown in Figure 10. 

854 5.1. Net Arctic Ocean Heat Output 

855 [ 47 ] Figure 10c shows an increasing trend in the total 

856 Arctic Ocean heating rate for the fall periods, while 

857 Figure lOd shows comparatively little change in the winter 

858 heating rate. Figures 10a and lOd show that for sea ice- 

859 covered regions, the net heating rate did not change mark- 

860 edly compared to the full Arctic Ocean domain in both the 

861 fall and winter. The heating rate over nonice-covered areas 


changed most dramatically in 2007 due to the larger amount 862 
of open water in that year (Figures 10b and lOe), increasing 863 
by nearly a factor of 5 from the previous years. Though ice- 864 
covered areas made up the dominant portion of the Arctic 865 
Ocean, the total heating rates were nearly equal over ice- 866 
covered and nonice-covered areas for the fall periods (with 867 
the exception of 2007). In 2004, 2005, and 2006 the net 868 
heating rate increased by 44%, 17%, and 12% from 2003, 869 
respectively. While in 2007 the large increase in nonice- 870 
covered areal coverage caused the total heating rate for the 871 
Arctic Ocean to increase by 300% from that in 2003. With 872 
the exception of the much later MA07 measurement time 873 
period, there was much less change in the winter time 874 
heating rates with a maximum change of 16% observed. 875 
[48] The results show an overall increase in the amount 876 
of ocean-atmosphere heat transfer in the fall periods. 877 
Section 4.2.2 showed that independent of changes in mete- 878 
orological conditions, thinning of the sea ice cover is 879 
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Figure 9. Winter time period ice thickness distributions, mean basal ice growth rates, and mean surface 
air temperatures for the'me-covercd east and west Arctic regions. 


880 responsible for up to a 40% increase in the net heat output 

881 in the ice-covered Arctic Ocean. However, this increase is 

882 small compared to the effect caused by changes in the ice 

883 areal coverage. The anomalously low areal coverage of sea 

884 ice in 2007 marked a turning point where the net Arctic 

885 Ocean heating rate became dominantly determined by the 

886 amount of ice-free area. 

887 5.2. Net Arctic Ocean Ice Production 

888 [ 49 ] The observed changes in sea ice thickness and ocean- 

889 atmosphere heat flux also lead to changes in the ice growth 

890 rate. Of particular interest is whether the observed losses in 

891 sea ice thickness and areal coverage led to a higher rate of 

892 ice production which could aid in the recovery of sea ice 

893 thickness and volume. 

894 [ 50 ] For sea ice-covered regions, the mean basal ice 

895 growth rates are shown in Table 6. Though basal ice growth 

896 varied with time depending on the surface air temperature 


and ice thickness distribution in a similar manner as the heat 897 
flux, Table 6 shows that a higher growth rate in the fall was 898 
generally followed by a lower growth rate in the winter and 899 
vice versa. The observed decreases in ice thickness may be 900 
due to a longer melt season as observed by Markus et al. 901 
[ 2009 ], increased oceanic heat flux as observed for the 902 
western Arctic by Woodgate et al. [ 2010 ], and/or increased 903 
ice export rather than due to changes in ice growth. These 904 
observations show that an expected increased basal ice 905 
growth rate associated with decreasing ice thickness did not 906 
largely occur over the 2003-2008 time period mainly due to 907 
associated changes in the surface air temperature. 908 

[ 51 ] The rate of ice volume production for ice-covered 909 
and nonice-covered areas is shown in Figure 10 , the pro- 910 
duction of ice can be seen to vary considerably from year to 911 
year. For the fall season ice-covered portion of the Arctic 912 
Ocean, the production of ice peaked in 2005 and 2006 due 913 
in part to the thinning of the ice cover and associated 914 
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Figure 10. Net ocean-atmosphere heating rate and ice volume production for the (a) ice-covered, 
(b) nonice-covered, and (c and d) total Arctic Ocean, (e) The dark colored bars represent the areal coverage 
of ice-covered regions, and the light colored bars represent the nonice-covered areal coverage. For the 
winter time periods, all regions are ice covered. The total area of the Arctic Ocean domain for all time 
periods in this study is 6.47 x 10 6 km 2 . 
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t6.1 Table 6. Basal Ice Growth Rate for Ice-Covered Regions During the Fall and Winter Seasons 3 




2003-2004 

2004-2005 

2005-2006 

2006-2007 

2007-2008 

t6.2 

Fall ice growth (cm month l ) 

10.1 (14.2) 

13.3 

18.1 

17.7 

12.4 

t6.3 

Winter ice growth (cm month -1 ) 

22.6 

19.2 

17.8 

14.7 

21.9 


t6.4 “The ON032 period is shown in parentheses. 


915 increased ocean-atmosphere heat flux discussed in section 4. 

916 In 2007, the production of ice in ice-covered regions 

917 reached the lowest point due to the high surface air tem- 

918 peratures and low ice areal coverage of the time period, 

919 while in nonice-covered areas the ice production increased 

920 by nearly a factor of 3 compared to the previous fall seasons. 

921 [ 52 ] For the full Arctic Ocean fall periods, the combination 

922 of ice production in ice-covered and nonice-covered areas led 

923 to a peak in the ice production in 2005 and a decrease in the 

924 following years. Despite the large increase in total ocean- 

925 atmosphere heat output in 2007, warm ocean and air tem- 

926 peratures kept the level of ice production near to that of 

927 2004. Thus, the 2007 ice minimum led to a greatly increased 

928 release of heat from the ocean to the atmosphere, however 

929 this increased heating rate did not lead to an increase in 

930 overall ice production because the ocean had yet to cool to 

931 the freezing point. The winter period ice production was 

932 much less variable, excluding the much later MA07 mea- 

933 surement period the ice production varied by less than 20% 

934 over the 2004-2008 time period. The winter time ice pro- 

935 duction variability was driven primarily by variability in the 

936 surface air temperature. 


937 6. Summary and Discussion 

938 [ 53 ] In this study we have combined ICESat freeboard 

939 retrievals with a snow depth model to estimate snow and sea 

940 ice thickness values for the Arctic Ocean during the 2003- 

941 2008 fall and winter time periods. The thickness data were 

942 used with meteorological data and a thermodynamic sea ice 

943 model to calculate the turbulent, radiative, and conductive 

944 heat fluxes, as well as the total ocean-atmosphere heat output 

945 and ice volume production for the Arctic Ocean. Sensitivities 

946 to the input parameters were determined and used to estimate 

947 the error in the calculated ocean-atmosphere heat fluxes and 

948 ice growth rates. The main factor affecting the uncertainty in 

949 our results was found to be uncertainties in the surface air 

950 temperature. Laser altimetry data was found to be particularly 

951 useful for determining the heat fluxes since the results are 

952 relatively insensitive to snow depth errors. 

953 [54] The heat flux and ice growth rates in ice-covered 

954 regions presented here are consistent with those from pre- 

955 vious observational studies conducted on multiyear ice. The 

956 advantage of the data sets used in this study is that they 

957 allow for estimates of heat flux over the entire Arctic basin. 

958 Also in agreement with the results of previous studies [e.g., 

959 Kwok et al, 2009; Giles et al., 2008; Maslanik et ah, 2007], 

960 this study shows that during the 2003-2008 time period the 

961 mean Arctic sea ice thickness decreased with much of the 

962 thickest ice (>3 m) being replaced by ice 0. 8-3.0 m thick. 

963 Variability in the calculated ocean-atmosphere heat flux and 

964 basal ice growth for ice-covered regions was primarily 

965 driven by changes in the surface air temperature as well as 

966 by the observed changes in the ice thickness distribution. 
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Heat fluxes during the fall periods were more sensitive to 967 
changes in the ice thickness distribution, with the eastern 968 
Arctic experiencing the greatest change in ice growth and 969 
heat flux due to changes in the ice thickness distribution. 970 
Taking variations in meteorological conditions into account, 971 
the fall period ocean-atmosphere heat fluxes were found to 972 
be greatly increased in 2005, 2006, and 2007 compared to 973 
2003 due to thinning of the sea ice cover. The winter time 974 
heat fluxes were much more impacted by changes in the 975 
surface air temperature rather than changes in the ice thick- 976 
ness distribution. Although the mean ice thickness decreased 977 
over the 2004-2008 winter time periods, the winter effective 978 
insulation did not largely change until 2008 at which time it 979 
experienced a large decline of nearly 1 m in effective sea ice 980 
thickness. The large decline in the winter 2008 effective 981 
insulation is also associated with an increase in the heat flux 982 
after differences in meteorological forcings are taken into 983 
account, though this increase is not as prominent as that 984 
observed in the fall and is within the estimated uncertainty. 985 

[ 55 ] For the whole of the Arctic Ocean, this study shows 986 
that increases in the net ocean-atmosphere heat output have 987 

urred due to thinning and area (volume) loss of the 988 
sea ice cover. However, a remaining question is: what 989 
Initude of changes to the surface air temperature have 990 
curred due to this decrease in sea ice volume and asso- 991 
ciated increase in the ocean-atmosphere heat flux? Surface 992 
air temperatures in the Arctic are highly variable so quan- 993 
tifying the impact of a changing sea ice cover on surface air 994 
temperatures is difficult [Serreze and Francis, 2006]. Serreze 995 
et al. [2009] show that decreases in the areal extent of Arctic 996 
are tied to increased surface air temperatures for the 1979- 997 
2007 fall seasons, but that this effect is not largely present 998 
during the winter season. The increased surface air tem- 999 
peratures in the fall were found to be due to a surface heating 1000 
source and attributed to an increased surface heat flux. This 1001 
study shows that over the 2003-2008 time period losses in 1002 
both ice thickness and areal coverage did indeed lead to an 1003 
overall increase in the surface heat flux. Despite large losses 1004 
in ice thickness and effective insulation, changes in ice areal 1005 
coverage were found to be the dominant factor in impacting 1006 
the surface heat flux. Most notably, the anomalously low 1007 
areal coverage of sea in the fall of 2007 led to an ocean- 1008 
atmosphere heat output nearly 3 times higher than that from 1009 
previous years. 1010 

[ 56 ] Serreze et al. [2009] also note that slight warming 1011 
may also be beginning to appear in the winter time. They state 1012 
this may be due to delays in autumn freezeup, but eventually 1013 
decreased ice extent and thickness in the winter will also 1014 
begin to play a role. Delays in autumn freezeup have been 1015 
observed by Markus et al. [2009]. However, this study shows 1016 
that though there was a decrease in the mean thickness and 1017 
amount of thick (>3 m) ice in the winter, these changes did 1018 
not lead to a large change in the ocean-atmosphere heating 1019 
rate since it is less sensitive to changes in the amount of 1020 
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1021 thick ice. It appears that a surface wanning signal associated 

1022 with a thinning sea ice cover could just be beginning to 

1023 emerge in the winter, but future observations will be required 

1024 to determine whether this effect becomes stronger and more 

1025 significant with time. 

1026 [ 57 ] Overall, these results show that the decreasing volume 

1027 of the Arctic sea ice cover has led to a decreasing ability to 

1028 insulate the atmosphere from the relatively warm underlying 

1029 ocean. This effect is currently most pronounced in the fall, 

1030 with the winter being less affected as the ice has sufficiently 

1031 thickened to a point where the ocean-atmosphere heat flux is 

1032 less sensitive to changes in the ice thickness. These increased 

1033 heat fluxes in the fall periods likely played a role in increasing 

1034 surface air temperatures in the Arctic. Though this data set 

1035 spans only 5 years, it was collected at a time when large losses 

1036 in sea ice thickness and areal extent were observed. The 

1037 continuation of large-scale sea ice thickness measurements 

1038 from future airborne and satellite missions such as NASA’s 

1039 Operation IceBridge and the planned ICESat-2 mission, as 

1040 well as ESA’s CryoSat-2 mission, will be vital to under- 

1041 standing future changes to the sea ice cover and its impact 

1042 on the climate. 

1043 [58] A major limitation in this study of the Arctic ocean- 

1044 atmosphere heat flux and ice growth rate is the irregular time 

1045 sampling and limited temporal availability of ICESat data. 

1046 Future satellite altimetry missions will maintain year-round 

1047 data collection for improved observation of year-to-year 

1048 variations. For the currently available ICESat data, it would 

1049 be useful to combine the observational data with model data 

1050 using an assimilation approach. Doing so would enable a 

1051 better understanding of reasons for the large losses in iqiij' 

1052 volume over the time period, how annual ice production was 

1053 affected by the observed changes, and how an increased 

1054 ocean-atmosphere heat flux from a reduced ice cover affected 

1055 surface air temperatures throughout the whole of the Arctic. 
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